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Приводится асимптотическое решение задачи гидродинамического расчёта радиального двухслойного пори- 
стого подшипника бесконечной длины с учётом анизотропии проницаемости пористых слоёв. Данное реше- 
ние основано на полных нелинейных уравнениях Навье — Стокса и уравнении Дарси. Рассматривается слу- 
чай, когда проницаемости пористых слоёв меняются непрерывно в радиальном направлении, а затем — слу- 
чай, когда оси меняются как в радиальном, так и в окружном направлениях. В результате найдено поле ско- 
ростей и давлений в смазочном и пористом слоях. Получены аналитические выражения для основных рабо- 
чих характеристик подшипника. Дана оценка влияния нелинейных факторов, а также анизотропии на основ- 
ные рабочие характеристики подшипника. Найдены условия, при которых радиальный неоднородный двух- 
слойных подшипник по несущей способности обладает свойством подшипника двойного действия. Установ- 
лены наиболее рациональные значения конструктивных, режимных и других функциональных параметров, 
определяющих работоспособность подшипника. 

Ключевые слова: пористый подшипник, двухслойная пористая втулка, несущая способность, коэффициент 
трения, коэффициент нагруженности. 


Введение. Как известно [1, 2], применение вкладышей из спечённых пористых металлокерами- 
ческих сплавов позволяет улучшить работу подшипников скольжения. Эти сплавы обладают ан- 
тифрикционными свойствами. Они работают с меньшим износом, чем изготовленные из цветных 
металлов, и позволяют создать режим жидкостного трения за счёт замены масла в порах. Гидро- 
динамический расчёт радиальных пористых однослойных подшипников — достаточно исследо- 
ванный вопрос. В работах [3—5] приводятся расчётные модели однослойных пористых подшип- 
ников бесконечной и конечной длины. Здесь проницаемость пористого слоя в радиальном 
направлении считается постоянной. В то же время на поверхности пористой плёнки, прилегаю- 
щей к смазочному слою, и в пористой втулке проницаемость задаётся скачкообразно. Это не поз- 
воляет написать непрерывные условия сопряжения при использовании многослойных пористых 
втулок, что является существенным недостатком приведённых в [3—5] расчётных моделей. До 
настоящего времени в известной нам литературе не представлено теоретическое исследование 
работ многослойных пористых подшипников с учётом анизотропии проницаемости пористых сло- 
ёв в радиальном направлении. Ниже приводится решение этой задачи для двухслойного пористо- 
го подшипника бесконечной длины с учётом анизотропии проницаемости пористых слоёв и нели- 
нейных факторов. 

Постановка задачи. Рассматривается установившееся течение вязкой несжимаемой жидкости в 
зазоре двухслойного пористого подшипника бесконечной длины с учётом анизотропии проницае- 
мости пористых слоёв в радиальном направлении. Предполагается, что подшипник неподвижен, а 
шип вращается с угловой скоростью ©. 

Поместим начало полярной системы координат в центре подшипника. 


* Работа выполнена в рамках инициативной НИР. 
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Тогда уравнения контуров шипа и подшипника (рис. 1) можно записать в виде 
1 : 
с :г=а(1+Н),г=Ьг=Ь+В,Н = 20050 - ЗЕ’, === <1 (1) 


где й — толщина двухслойного пористого слоя, а — радиус шипа, 6 — радиус подшипника, е — 
эксцентриситет. 





Рис. 1. Схематическое изображение радиального двухслойного пористого подшипника 


Проницаемость пористых слоёв зададим таким образом, чтобы на границе раздела слоёв 
они совпадали и кроме того удовлетворяли уравнению Дарси для каждого слоя. 
Пусть 


ми № 
К=Ае °^, К=Ае °^. (2) 
Здесь А — заданная постоянная величина, характеризующая проницаемость на границе раздела 
слоёв; Ли, Л› — характеризуют распределение проницаемости слоёв в радиальном направлении. 
Основные уравнения и граничные условия. Будем исходить из безразмерных уравнений 
движения вязкой несжимаемой жидкости, уравнения неразрывности и уравнения Дарси [2]: 





























ди, и, д, и? 1 ор 0% 10, 10% 06 Ш 29, 
Кеч, —^ + +9, = - + А 5 5 й 
ог г 09 Г (1-а) 0г 0’ Гог г“ 09 02° г’ Г’ 00 
2 2 2 2 
Ве о Сы: у.е | _ 1 о В Е 
г г 00 Г (1-а)' г 99: >В“ РОО 5800" ОЕ” СР СР -00 (3) 
ди, 10% ц, 0Ф 10Ф 10Ф Л 
+ +—^=0, = р 
де г 00 г д” ог г 08 г д! 
0, 10% ц, РЕ 10Е 10+ МЛ, 0Е 
+ +—=0, + + = | 
0г г 00 г д"? гог г? 06? г дг 


Здесь ч,, и, — безразмерные компоненты вектора скорости; р,Ф, Е — гидродинамическое дав- 
ление соответственно в смазочных и пористых слоях. 
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Размерные величины г, ч,, ц%, р, Ф,Е, К, связаны с безразмерными Г, 0,,0,, р, Ф, Е, К, 


соотношениями: 
г=ЬГ, ч, =Оад,, цч=Фа0,, р = рае: Б, 
(6-а) 





(4) 
ф- ы@ар ФЕ - ы@аБ = 
Ь-а Ь-а 
В дальнейшем знак ^ у безразмерных переменных опускается. Рассмотрим граничные 
условия, при которых решается система уравнений (3). 
Пористая поверхность, прилегающая к смазочному слою, характеризуется следующим 
образом: 


1. На контуре с, выполняется условие непрерывности гидродинамического давления, а 


Е, К! = АК,, 1=1,2. 





компоненты вектора скорости определяются законом Дарси: 
К, ОФ 


т ы ог 





г=Ь 

2. На поверхности пористого вкладыша, прилегающего к непроницаемому корпусу, 
ОФ [0г=0. 

3. Выполняется условие замкнутости смазочного слоя (0) = р(2п), Ф.(г, 0) = Ф\(к 2п), 
Ф2(г, 0) = Ф2(+, 2п). 

4. На поверхности шипа выполняются следующие условия 





у, (а+аН)=ц,|_ +5 аН +...=-=510, 
г=а аг и 


чз (а+аН) = |, +5 аН+...=1. 


ог 
5. На границе раздела слоёв — равенство давлений, а также условие 
0Ф = и ФЕР. (5) 
ог ОГ 


Асимптотическое решение задачи. Решение системы (3), удовлетворяющее выше приведён- 
ным граничным условиям, будем искать в виде 


р=У =", у, Уи =“, че У", Ф-УФ, =", Е УЕ =“ | (6) 
Ю Ю Ю Ю Ю 


Для определения коэффициентов разложений (6) с точностью до членов О(=^) придём к следую- 
щей системе уравнений и граничных условий к ним: 


ве ди, * 1 бр 9% 1 Ш ве 44 а |. 0 19% 6 



































+ Й 
а т (1-а) 4 4 Гч “а: Е ар” ва. 27 @) 
Чи Ш _ о а?Ф, ‚ 1аФ, __ № а®, а. ‚14 __ № 
г “а га О а га 
2 2 
Ве| и, ди, +и, аи, - № ди, 2, и т | др, й д и Хх 1 ди, р . д о — - ди, | 
ог фр г 00 Г (1—а) ог ого Рог г 
2 
Ве би А аз 10 , ти Ш, _2 ди ‚ (8) 
[в/в 00 Г (1-а)' Г 0 * Гог 1? 00 г? 00' 
ди 10 0°Ф, 10, 1 0Ф, М0Ф ФЕ 10Е ТОЕ № 0Е 
+ +— =0, ы—% 5 о о : 
ог г 00 г ог г ог и 00 г г д’ го г’ 00 га! 
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Системы уравнений (7) и (8) решаются при следующих граничных условиях: 
и (а) =0, и, (а) =1, м (1) =-+(1-а).Ф, (1), 2, (8, =Ф. (В,), 





9 
Е, (В, ) =Ф, (В), ч, (1)=0, Ф, (1) = р. (1), Е'(В,)=0. (9) 
и, (а, 0) =-эт6, ц, (а, 0) =-— (а)асоз6, и, (,6) = 41 -в) 57 ро, (,6) =-ч4(1-в)%5 т 
В (6) -Ф(,6), РВ), (7,0) = (7,21), <, 8.8) (8,6), 9..6) АВ). 
Здесь В, =1+ Е 85 -1+р ф= а ка= г 


Решение системы (7), удовлетворяющее граничным условиям (9), легко находится непо- 
средственным интегрированием. В результате получим 








и =0, р, = Ве (1 с) [* а +с, Ф, = сопз = р, (1), 4 = ан =), Е, =соп8. (11) 
Исходя из вида граничных условий, решение системы (8) для первого приближения будем искать 
в виде 

и =и, (г) 059+, (г) тб, и, =ц,, (г) с0$0 +ц,, (г)зт6, р, = р. (г)со$6 + р,, (г) те, (12) 
Ф, =Ф, (г) с050+Ф,, (г) т, Е =А, (г)со$0 +2, (г) 56. 


Подставляя выражение (12) в уравнения (8) и граничные условия (10) и приравнивая выражения 
слева и справа при со50 и т, получим следующие уравнения и граничные условия: 





















































ие Ти И а й р’. +Ве| ши’ +и'и и а 
11 Г 11 [2 11 г? 12 (1 -@) 11 0-11 011 Г 12 Г 11 |7 
ВИ 2 1 В | р И 2%, 
ты. (-а} ВА | 
из и, = : —. + Ве или’, + и’, + би, + и + Фи 
11 г 11 2`и г? 12 (1 -а) г 12 0-11 0-11 г 12 г 11 г 11 |7 
1 2 2 1 1 и и у 
+20 и и р. +Ве ши. +иъи +, +-И (13) 
12 г 12 р 12 22 11 (1 -а) г 11 0-12 12-0 г 11 г 12 г 12 |! 
г. и а и 
И ++ =0, Ш --и,+-В=0, 
Ге г г Г 
" 1 " 1 Л, ’ п" 1 п" 1 Л, ’ 
Ф!, + г Ф!, 2 Ф,, г >, Ф,, + г Ф,, г? Ф,, а г Ф,, 
1 1 Л 1 1 
Е" 1 Е" Е 2 Е’ Е" 4 Е" Е 2 Е’ 
11 г 11 22 11 г 117 12 г 12 г? 12 г 12 


и (а) =0, и, (а) =-1, ц, (а) = м, (а)-а, ч,, (а) =0, и, (1) =-ф(1-а)-Ф!, (1), 
и (1) =-4+(1-а).Ф,, (1), ц,, (1) =-ф(1-а).Ф,, (1), ч,, (1) =%(1-а).Ф,, (1), 
Р, (1) — Ф,, (1), Р, (1) ы Ф,, (1), Е. (В») =0, Ра (В,) =0, Ф,, (В,) = ЕР [В 
Ф, (В,) =А. (В,), Ф (В,) = Ал (В), Ф, (В,) = (В,). 
Заменяя в выражениях системы (13) и (14) производные слагаемые конечноразностными пред- 


ставлениями, получим систему алгебраических уравнений, которая решается методом Гаусса — 
Зейделя. 


(14) 
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Определение основных рабочих характеристик подшипника. Определив поле скоростей и 
давлений в смазочном слое, можно перейти к определению основных рабочих характеристик 
подшипника. Для составляющих вектора поддерживающей силы А, и А), а также для момента 
трения получим следующие выражения: 

=пра*Бо р =п. ра? Е 2азрОп 
= > р, (а), А, = > р. (а), М, = Ь 
(6-а) (6-а) 


Коэффициент нагруженности с, коэффициент сопротивления вращению & и коэффициент трения Ё 
являются основными рабочими характеристиками рассматриваемого подшипника. Они определя- 
ются по формулам: 


+ (а) +0 (=). (15) 


Хх 








о ее Е (16) 


21аро 2180 '’1-а с 
Кроме того, на работоспособность рассматриваемого пористого подшипника влияют следующие 
параметры: 


.. 
№ п—1 


— постоянная проницаемость стенки вкладыша К“ = Ае И. характеризуемая параметром 
ф = КР(Ь-а); 

— радиальный относительный эксцентриситет = = Е/(Ь-а); 

— толщина вкладыша, характеризуемая безразмерной величиной В = (6+1)/В; 

— отношение толщины пористых слоёв /ь/№!; 

— длина подшипника / 

— число Рейнольдса (Ве); 

— параметры ^Л, и ^, , характеризующие распределение проницаемости в радиальном направле- 





нии в подшипнике. 
Следует отметить, что в случае учёта анизотропии проницаемостей пористых слоёв в ра- 
диальном и окружном направлениях имеем 


ми ит 
К =Ас0$00е ^*, К, =Асо$0е ^”*. 


Осредняя по окружности выражения для К, и К, , получим 


ми № 
© А чинов ‚К -_А чп2пюе* "8. (17) 
т 216 * 21° 


Полученные выше аналитические выражения для основных рабочих характеристик подшипника, 
определяемые формулами (15), в рассматриваемом случае остаются без изменения, если здесь 





Ь 

МАт—— . Ь 
АБе` 2*® Абзт21м Мп 
выражение для у =-——— заменить выражением ц = ^^^ е т 


(6-а) 210 (6 -а)} 
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Численный анализ полученных результатов. 








ФЗ 41 В 














13 
Рис. 2. Зависимость безразмерной несущей способности Рис. 3. Зависимость безразмерной несущей способности 
подшипника от параметров В и @ (при значениях подшипника от параметров В и @ (при значениях 
параметра @>1), Л! = Л» = 0, а = 0,998: параметра &<1), Л! = Л» = 0: 1 — без учёта сил инерции; 
1 — без учёта сил инерции; 2 — с учётом сил инерции 2 — сучётом сил инерции 
$ 27 
1-а 
0,7 
0,6 
0,5 
0,4 ыы 
0,3 
0,2 10 
0,1 
103 102 107 1 ы 103 107 107 1 у 
Рис. 4. Зависимость коэффициента нагруженности Рис. 5. Зависимость коэффициента трения от параметра 4: 
от параметра ц, а = 0,998: 1 — А >1, п <1, ® = 0; 1— лы >1, В, <11 ®=0;2— Л, = № =0, ®=0; 
Л, Р, Л, В, 
А В Л В 
2— № =№=0, = 0; 3 — + <1, -2>1, 9=0 3— 1 <1, -2>1, Ф=О 
А х р йа” 


2 1 2 


Выводы. 
1. Поддерживающая сила не перпендикулярна линии центров при: 


— Л, =Л, =0, т. е. проницаемость однослойного пористого вкладыша определена формула- 
ми (17) (К, = К, =0) 


— учтены силы инерций. 
Значение К, — составляющей поддерживающей силы на 5—10 % больше, чем при Ке = 0. 





2. Если Л, =Л, =0, К, =К, =0, При & - подшипник по несущей способности обладает 


свойством подшипника двойного действия. При &<1 значение несущей способности стабили- 


зируется. 
3. Если Л, = Л, , с увеличением значения параметра ф значение коэффициента нагружен- 


ности уменьшается. Наиболее резкое снижение значения коэффициента нагруженности наблюда- 
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Л в 

ется при ш=10?. При о я >1 значение коэффициента нагруженности больше, чем при 
2 1 

Л в 


—>1, -2*<1, 
^, 1 


4. Если Л, =, ‚ с увеличением значения параметра Ш значение коэффициента трения 
снижается. Наиболее резкое снижение значения коэффициента трения наблюдается при ш=10?. 
При значениях у =10" зависимость коэффициента трения от параметра ф минимальна. 


Таким образом, учитывая анизотропию проницаемости пористых слоёв, определяемых 
формулами (17), можно обеспечить повышенную несущую способность подшипника при доста- 
точно низком значении коэффициента трения. 
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НУРВОРУМАМТС САЕСИЕАТТОМ ОЕ ТМ/О-ЕАУЕВ РОВОЧ$ ВЕАВТМС ОЕ ТМЕТМШТЕ 
ГЕМСТН М/ТТН АССОЦМТ РОВ РЕВМЕАВИЛТУ АМТЗОТВОРУ ОЕ РОВОЦ$ 1АУЕВ. 
АМО ТМЕВТТА РОВСЕ$” 


К. $. АКВуегоТуеу, М. А. МикКцаад?е, \. $. Моудогодома, Т. $. СпегкКазоуа 
(Козоу Зае Тгапзро": Упмег®Ку) 


Тре азутрюйс зо/ийоп {0 а рго Мет о! {те пуагодупатис са/сшайоп о! пе габ/а! {ио-Шуег рогоиз Беата о! тЕ- 
пе 1епдёй и! ассоипЕ Гог {пе реттеа Шу ап/зо{гору о! рогоиз !ауегс [5 ргезещеа. ТПе зо/иНоп 15 Базе оп и! 
поптеаг Мамег — ЗюКез едиавоп5 апа Рагсу5 едиайоп. ТВе сазе ипеп {те рогоиз !ауег регтеа Шу спапдез 
сопйпиоиз/у т {те габа! @гесйоп, ап еп — {те сазе итеп ахез свапде Бой т габа! апа тп сгситГегепва! 
!гесНоп5 аге сопз/СЕгед. Аз а гезий, а уе/осйу апа ргеззиге Неа т {те ибтсапЕ апа рогоиз [вуегс 15 Ююипа. Апа- 
/уйса! ехргезюоп5 Юг пе тат регогтапсе (а о! {те Беатпд аге обатеа. ТВе ейесЕ о! поп/теаг !асюг$ апа 
ап/зо{тору оп {те Беаппд Бах регогтапсе /5 а55е55е4. Тре сопа!юоп5 ипаег иР/сй а габа/ поп-ипйогт &ио- 
[ауег Беата ассогатд {0 {те саггутд сарасйу отегх пе дои Ме асвоп ргорепу аге дЕегттед. Тве тозЕ гаНопа! 
уа/иез о! пе дез!одп, орегайпа, апа о{Тег Гипсвопа! да дейп/та {пе Беагта орега Шу аге езаБ/зред, 
Кеуигога&: рогоиз Беата, В/о-1уег рогоиз р/ид, Беата сарасйу, Е7своп соейсепЕ !оаатд соейсепЕ, 


* ТНе гееагс {5 аопе мт {Пе Нате оЁ {Не паерепдепе В&р. 
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